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We investigate an ideal gas in a time-dependent external trapping potential. We use the 
Boltzmann equation with the relaxation time ansatz to explore the time-dependent energy of an 
adiabatically isolated system. In particular we are interested on the dissipation during a potential 
change along a given protocol with finite velocity. The role of the relaxation less evolution as a 
limiting case is studied: starting from an equilibrium distribution and for small times the energy 
of the gas with relaxation is always smaller than that of the relaxation-less. This means that 
relaxation processes show an ambivalent behavior: on the one hand entropy production, but on the 
other hand reduction of dissipation by driving back the system into equilibrium. 
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I. INTRODUCTION 


Dissipation of energy in gases or liquids due to a mov¬ 
ing body is at the core of physical theories such as hy¬ 
drodynamics, statistical physics or kinetic theory. Some 
interesting new results at a fundamental level, regard¬ 
ing the connection of equilibrium and non-equilibrium 
free energy differences, as the Jarzynski equation, or the 
relation between dissipation and the asymmetry of time- 
reversal processes are known[Il-Q. These results have 
been extended from classical to quantum mechanical sys¬ 
tems In this work we investigate the dissipation 

process in gases using the Boltzmann equation (BE) with 
a time-dependent external potential. The main moti¬ 
vation for considering this problem was the ambivalent 
behavior of elastic collisions during the potential defor¬ 
mation; on the one hand collisions do not change the 
energy of the colliding particles, and this is also true for 
a time-dependent potential, since the collisions are de¬ 
scribed as pointlike in phase space. On the other hand it 
is immediately clear that the energy of the gas in a time- 
dependent potential must depend on the collisions during 
the potential change. It makes a difference in the distri¬ 
bution function if in the BE collisions are present or not. 
The first impulse to investigate such problems using ki¬ 
netic theory was the study of Bose-Einstein condensation 
in trapped gases, where one can reach the condensation 
of a Bose gas by changing the trapping potential [1, [l^ . 
The paper is organized as follows: Sec. HIl gives a short 
overview of the theoretical elements and formulas used to 
derive the main results represented in Sec. nni Sec. nv] 
gives a numerical example and Sec.|V]conclusion and out¬ 
look. The appendix contains a proof of the H-Theorem 
for time-dependent potentials and the relaxation time 
ansatz in the BE. 
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II. GENERAL APPROACH: KINETIC THEORY 

First we present some general equations of kinetic the¬ 
ory used in this paper. The central quantity to de¬ 
scribe the system is the time-dependent one particle 
distribution function /(p,x, t). The distribution func¬ 
tion depends on the one particle phase space coordi¬ 
nates {p, x} (momentum and space coordinates, some¬ 
times suppressed to simplify notation) and time t (also 
suppressed sometimes). / is normalized to the number 
of particles N, 


J f{p,K,t)dT = N . (1) 

Here dT = d^a;d^p/(27r)^ denotes the volume element in 
the 6 dimensional phase space (we omit always Planck’s 
constant h and the Boltzmann constant ks)- 


A. Kinetic equation 


We consider a gas in a time-dependent external po¬ 
tential U{'x.,t). The time evolution of the distribution 
function is given by the Boltzmann equation (see for ex¬ 
ample [ni) 

f + (2) 

Here {df /dt)coii denotes the collision integral, which we 
approximate by the relaxation time ansatz. 


f-fa 
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( 3 ) 


Because the external potential is time-dependent 
the equilibrium distribution /g(p,x, t) becomes time- 
dependent as well. 


fg{p,x,t) = exp{-l3{t)[H{p,x,t) - fj.{t)]} . (4) 
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Here ^{t) is the time-dependent chemical potential, 
/3{t) = l/T{t) corresponds to the time-dependent inverse 
temperature T{t) and 

= ^ + U(x,t) , ( 5 ) 

is the one particle Hamiltonian for particles with mass 
m and without inner degrees of freedom. The equilib¬ 
rium distribution Eq. o does not contain a macroscopic 
momentum or angular momentum because we restrict 
ourselves to spherical symmetric potentials to make the 
discussion transparent. The main result Eq. m is not 
affected by this restriction, as the derivations show, since 
the explicit structure of fg does not enter. So let us as¬ 
sume that the number of particles and the energy for a 
fixed time are the only conserved quantities, which al¬ 
low the determination of /3(t) and ^{t) in Eq. Q. In 
particular we have the relation 

y/ 9 (P,x,t)dr = y/(p,x,f)dr , (6) 

ensuring conservation of the number of particles. The 
conservation of energy E(t) for a fixed time is expressed 
as 

^^(0 = y -ff(P,x,f)/g(p,x,t)dr = 

y iJ(p,x,f)/(p,x,t)dr . (7) 

Eq. m with an initial condition (in our case a Maxwell- 
Boltzmann distribution) and the two conditions Eq. (|5|) 
and Eq. 0 are the basis for our analytical and numerical 
studies in Sec. IlYl We remark that Eq. ([5]) due to Eq. (|5|) 
and Eq. o becomes non linear, in contrast to the Focker- 
Planck equation, which is a linear evolution equation. 

B. Entropy 

A central quantity is the nonequilibrium entropy of the 
gas. For a given one particle distribution function it is 
computed as m 

S{t) = - J /(p, X, t) ln[/(p, X, t)]dr . (8) 

One can proof the validity of the H-Theorem, S{t) = 
— J \n{f){df/dt)coiidT > 0 with the relaxation time 
ansatz and a time-dependent external potential also. For 
details see the appendix App. Here we need an addi¬ 
tional expression for S{t) for gases not far from equilib¬ 
rium. We split the distribution function as / = /g -f Sfg, 
where Sfg denotes the small deviation of / from the equi¬ 
librium. Expansion of the logarithm in Eq. ([8]) up to 
second order and some lines of algebraic simplifications 
yields 

S{t) = S,g{t) - I ^dT . (9) 


Here Seq{t) = f5{t)E{t) — l3{t)^{t)N is the time- 
dependent equilibrium entropy, reached when Sfg = 0 
[1^ . Due to relaxation Sfg{t) becomes smaller and 
smaller and at the end (for a fixed potential) S{t:^T) = 
Seq- 

III. DISTRIBUTION FUNCTION: 

PERTURBATION THEORY 

To proceed we have to compute the distribution func¬ 
tion. Because the analytical integration of the kinetic 
equation for time-dependent potentials (even with the 
simple relaxation time ansatz) in general is not possible 
we look for a perturbation theory approach. There are 
two interesting limiting cases as starting points for this 
attempt: an infinitely slow potential change, correspond¬ 
ing to the adiabatic time evolution of the distribution 
function, and the relaxation less time evolution [t ^ oo 
in Eq. 0], corresponding to the Liouville evolution (de¬ 
noted as /l), respectively. We choose here the second 
way, where the system has a large relaxation time, com¬ 
pared with the time t, for which we are interested, so 
that t/r ^ 1. An analytical integration of the time- 
dependent collision less BE is still impossible in general, 
but here we do not need the solution explicitly. 

A. First order expression 

The lowest (zeroth) order of such an expansion is the 
time-dependent Liouville function f^. We split the true 
distribution function into two parts, 

f = fL+SfL (10) 

where /z,(p,x, f) is the solution of the collision less BE 
and SfL is the deviation from the actual distribution. If 
we enter with this ansatz into Eq. 0, we get for the first 
order correction 

^/l’(P, x,t) =-- f dt'[fL{p,x,t') - f^°'>{p,x,t')] . 

T Jo 

( 11 ) 

This first order expression is sufficient as long as the ratio 
t/r is small. The equilibrium distribution fg^'^ is com¬ 
puted using //, on the right hand side of the conditions 
Eq. (jni), Eq. ([T]). We do not require the computation of 

but we know that it does not depend on the relax¬ 
ation time T. 


B. Energy 

The distribution function Eq. m with Eq. (HU allows 
us to clarify the influence of r on the dissipated energy. 
We consider the begin of a potential change, where the 
distribution function deviates weakly from equilibrium. 
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The starting point is Eq. where we replace S (t), using 

Eq. (IA5I) , by 


S{t) = s{0) + 1 f ^dr. (12) 

^Jg 

Now we use / from Eq. where we approximate the 
time integral in Eq. (EH) by a trapezoidal rule. After 
some simplifications we obtain up to the order (t/r)^, 


S,git)=S{0)+ 1 



Ul - 


2/. 


( 0 ) 


dr . (13) 


The dissipated energy (heat) is given by T{t)[Seq{t) — 
5'(0)]. Eor small times the choice T{t) = Tad{t) is rea¬ 
sonable and we get [13 


E{t) — Ead{t) +Tad{t) 1 





2/^ 


( 0 ) 


dr. 


(14) 


For T —>■ oo we obtain from this expression 


EL{t) := J HUdV = Eadit)+Tad{t) 

and further (for t/r <C 1) 


(/L-/f )^ 

2/i°^ 


dr. 
(15) 


E{t) < EL{t) . (16) 


This result shows that the time-dependent energy El 
represents an upper bound of the energy of the gas, in¬ 
dependent of the characteristics of the potential. Re¬ 
laxation favors minimization of dissipation. This is re¬ 
markable because the Liouville evolution itself conserves 
entropy S{t) = — J Jl ln(/i)dr = S{0). For larger times 
Eq. (fT6)l can be violated, however for slow potential 
changes it must fulfilled, because the true distribution 
functions / and fg approaches the adiabatic distribution 
function fad- Because the dissipated energy is given by 
Ediss — E Ead One obtains Ediss — El Ead- 


IV. NUMERICAL EXAMPLE 

For a model potential we consider a one dimensional 
harmonic oscillator with time-dependent frequency, as 
widely used to investigate dissipation of Brownian par¬ 
ticles and the Jarzynski equation [i, [13. An in¬ 
crease/decrease of the frequency leads to a compres¬ 
sion/expansion of the gas. We integrate numerically the 
BE with the relaxation time ansatz and monitor the en¬ 
ergy (and entropy) of the gas during the potential change. 
The numerical method is based on a discretization of the 
phase space into equal sized cells and a numerical inte¬ 
gration of the resulting ordinary differential equations. 
To obtain the time-dependent equilibrium distribution 
fg, we compute for each time step the integrals Eq. 


(a) Energy 




FIG. 1: (color online) (a) Energy of the gas during the com¬ 
pression and expansion. The frequency ratio is Aui/ujq = 0.5. 
The green line corresponds to the Liouville evolution (El) 
and the blue line to a evolution with finite relaxation time 
{r/tf — 0.5). The red dotted line is the adiabatic limit and 
in black-dotted it is shown the result from the perturbation 
theory, (b) shows the dissipated energy for a time regime 
where Eq. EH is valid. In (c) we have plotted the entropy of 
the gas and the approximation according Eq. EH- For small 
times El is the upper bound for the energy. 


and Eq. m and solve numerically the resulting equations 
for P and /r (where the values for the quasistatic poten¬ 
tial change are used as a starting guess). This procedure 
runs well and there are no numerical problems which we 
should report. To be sure that energy conservation (not 
only number of particles) for a fixed potential is fulfilled, 
the used numerical limits of the phase space must be 
large enough. The external potential has the form 

U{x,t) = , uj{t) = ujQ + Aujsm{2Trt/tf) . 

(17) 

Here tf characterize the switching time and likewise de¬ 
notes the final integration time. We set t/ = 1 as our 
time scale. We choose wot/ = 1 and jAw/wol = 0.5, so 
we are far away from adiabatic conditions. The mean free 
fly time of a particle is on the order ^ 1/{u}qN), where we 
have choosen N = 1250 (and the ratio initial tempera¬ 
ture to initial frequency Tq/ojq = 200). For Bose gases at 
low temperatures and densities as H and Na in magneto¬ 
optical traps the relaxation time differs dramatically due 
to the different s-wave scattering cross sections (by a fac¬ 
tor of ^ 10^l[l3- For an adiabatic potential change the 
temperature of the gas is given by Tad/T\i = w/wq and 
the energy Ead/Eo = w/wq, resoectivelyfl^. In Fig. [T] 
we show the energy/entropy for one compression - expan¬ 
sion chicle as a function of time for a gas with specified 
relaxation time. The lines for the energy reflect the time- 
dependent frequency, whereas the entropy increase as re¬ 
quired from the H-theorem. We have also plotted the 
perturbation theory expressions Eq. (ICT) and Eq. (fT^ to 
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Energy 



FIG. 2: (color online) Energy of the gas during the expansion 
and compression cycle, Acj/ojo = —0.5. Although it is difficult 
to see the energy corresponding to the Liouville evolution 
for small times is always greater than all others, in accordance 
to Eq. (IT6|l . 



FIG. 3: (color online) The same potential deformation as in 
Fig. m but 3.5 times slower as there (with r = tf). In (a) it 
is shown the initial Maxwell-Boltzmann distribu tion in phase 
space [in units {P,X) = (p/v^mZTo, a;y^mwJ/Tb)], (b) shows 
the partially relaxed distribution a.t t = tf and (c) shows the 
Liouville distribution /i(t/). In (d) we have plotted again 
the energy of the gas (blue line) to demonstrate the validity 
of Eq. (1161) for the entire time. 


demonstrate the result Eq. (fTfij). Fig. [2] demonstrates the 
same behavior of the gas for the opposite chicle, where 
first the frequency uj decrease and hereafter increase (ex¬ 
pansion and compression). Here we have plotted the 
time-dependent energy of two gases with different relax¬ 
ation times and additional Eadit) and For small 

times Ei,{t) is an upper bound of the energy. Fig.|3]shows 


the distribution function, obtained by the numerical inte¬ 
gration. The initial Maxwell-Boltzmanm distribution (a) 
and the collision less evolved function (c), as well as an 
example with finite relaxation time (b) the are plotted. 
In this case (slower potential change as in the previous 
example) the energy of the r 7 ^ 0 system is smaller as El 
for the entire switching time. 


V. CONCLUSION AND OUTLOOK 

In this paper we have discussed energy dissipation in 
gases in time-dependent trapping potentials. Using the 
Boltzmann equation with the relaxation time ansatz we 
have shown that the dissipated energy depends on the 
relaxation time. In particular we have demonstrated us¬ 
ing perturbation theory, that for small times (compared 
with the relaxation time) the relaxation less time evolu¬ 
tion always leads to a higher energy than that with finite 
relaxation time. On this time scales relaxation decreases 
dissipation. The derivation shows that this statement is 
fairly general: it does not depend on the shape of the 
external potential. It would be interesting to make some 
numerical studies with the full collision integral of the BE 
to check if the results of this paper are still valid dill. 


Appendix A: H—Theorem for time—dependent 
potentials and the relaxation time ansatz 


The time derivative of the entropy in the relaxation 
time ansatz is equal 

Sit) = J Hf)^^dT . (Al) 


If we set inside the logarithm f = fg + 6fg we obtain 

J-h. 


J lnif)^-^dT = i y /ln(/)dr - 


Iff, In 




dr. 


(A2) 


Now we split the argument of the logarithm in to two 
parts. 


Sit) = 


Segit)-Sit) 1 


--//.in 


1 -b 




dr , (A3) 


where we have defined the equilibrium entropy Seqit) = 
— f /gln(/g)dr. Because fg is the equilibrium distribu¬ 
tion, corresponding to the same number of particles and 
energy as /, it follows Seqit) — Sit) > 0. The last inte¬ 
gral in Eq. (IA3I) is positive for all t. If a; > 0 is a real 
number, then ln(a:) < a; — 1 is always true. Using this for 
the integral above we obtain (/g > 0 ), 


/gin 


1 + 14 

Jg J 


dr < y SfgdT = 0 


(A4) 
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With the — sign in front of the integral in Eq. (IA3I) we 
conclude that S >0. 

In the text we use an expression for S for small Sfg. One 
can obtain this expression by insertion oi f = fg + Sfg 
into the logarithm in Eq. and a series expansion, 

S{t) ^- [ ^dr . (A5) 

T J fg 
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